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Abstract 

We present in closed form some special travelling-wave solutions (on the real 
line or on the circle) of a perturbed sine-Gordon equation. The perturbation of the 
equation consists of a constant forcing term 7 and a linear dissipative term, and 
the equation is used to describe the Josephson effect in the theory of superconduc- 
tors and other remarkable physical phenomena. We determine all travelling-wave 
solutions with unit velocity (in dimensionless units). For \j\ < 1 we find fami- 
lies of solutions that are all (except the obvious constant one) manifestly unstable, 
whereas for \j\ > 1 we find families of stable solutions describing each an array of 
evenly spaced kinks. 
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1 Introduction and preliminaries 



The scope of this communication is the determination in closed form of some special 
solutions of the class of partial differential equations 

<Ptt ~ Pxx + sin ip + aupt + 7 = x £R, (1) 

parametrized by constants a > 0, 7 G R, more precisely the determination of the 
travelling-wave solutions cp(x,t) = g(x — vt) with velocity v = ±1. 

This equation (here written in dimensionless units) has been used to describe 
with a good approximation a number of interesting physical phenomena, notably 
Josephson effect in the theory of superconductors [6], which is at the base pQ of a 
large number of advanced developments both in fundamental research (e.g. macro- 
scopic effects of quantum physics, quantum computation) and in applications to 
electronic devices (see e.g. Chapters 3-6 in [2]), or more recently also the propa- 
gation of localized magnetohydrodynamic modes in plasma physics [9]. The last 
two terms are respectively a dissipative and a forcing one; the sine-Gordon equation 
(sGe) is obtained by setting them equal to zero. 

The sGe describes also the dynamics of the continuum limit of a sequence of 
neighbouring heavy pendula constrained to rotate around the same horizontal x- 
axis and coupled to each other through a torque spring [8] (see fig. [[]); (f(x,t) is the 
deviation angle from the lower vertical position at time t of the pendulum having 
position x. One can model also the dissipative term —apt of (JTJ) by immersing 
the pendula in a linearly viscous fluid, and the forcing term 7 by assuming that a 
uniform, constant torque distribution is applied to the pendula. This mechanical 
analog allows a qualitative comprehension of the main features of the solutions, e.g. 
of their instabilities. The constant solutions of (JTJ) exist only for I7I < 1 and are, 
mod 2ir, 

(f s (x, t) = — sin -1 7, <f u (x, t) = sin -1 j+ir. (2) 

If I7I < 1 the former is stable, the latter unstable, as they yield respectively local 
minima and maxima of the energy density 

h ._£ + !& + 1<p _ 00BVm (3) 

In the mechanical analog they respectively correspond to configurations with all 
pendula hanging down or standing up. If 7 = ±1 ip s = ip u = =F7r/2 (mod 2n), 
which is unstable because it is an inflection point for h. 

In [5] we have performed a detailed analysis of travelling- wave solutions of (pQ). 
We briefly recall the framework adopted there and some of the results. Without 
loss of generality we can and shall assume 7 > 0: if originally 7<0, we just need to 
replace ip — * —ip. Moreover, space or time translations transform any solution into 
a two-parameter family of solutions, so one can choose any of them as the family 
representative element; for travelling- wave solutions this reduces to translation of 
the only independent variable. In agreement with the conventions adopted in [5], 
we specify our travelling-wave Ansatz as follows: 

£:=±x-t, l p(x,t)=g(0-7r. (4) 
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Figure 1: Mechanical model for the sine-Gordon equation, (a) Spring, (b) solder, (c) 
brass, (d) tap and thread, (e) wire, (f) nail, (g) and (h) ball bearings, (i) base (After A. 
C. Scott [8], courtesy of A. Barone, see []]) 



Replacing the Ansatz in (TjQ) one obtains the first order ordinary differential equation 

ag = 7 — sin g. (5) 

We have already recalled the constant solutions. If g' is not identically zero, by 
integrating d£ = adg/^ — sing) one finds 




Co 90 



separately in each interval in which g' keeps its sign. This allows to determine the 
solution implicitly, namely the inverse £,(g). 

If 7 < 1, as g approaches respectively sin 7 or it — sin _1 7 (mod. 2ir) the de- 
nominator of the integrand goes to zero (linearly if 7 < 1, quadratically if 7 = 1) 
while keeping the same sign, and therefore the integral diverges, implying that the 
corresponding £ respectively goes either to ±00, or to =poo [5]. In either case the 
range of £(<?) is the whole R, implying that by taking the inverse one obtains g(£) 
already in all the domain. If 7> 1 the denominator of the integrand is positive for 
all s S R, so that the solution g is strictly monotonic and linear-periodic, i.e. the 
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sum of a linear and a periodic function, and 

f 2w ds 

<?(£ + H) = 5 (£) + 2vr, E := a / — . (6) 

Jo 7 - sin s 

Denoting as ip the corresponding solutions with £ :=±r — t, by © this implies 

^(x+E.t) = Cp ± {x,t)±2Ti. (7) 

This behaviour is illustrated in fig. [2] by a picture of the corresponding configuration 
for the mechanical model of fig. [TJ 

ip^ can be interpreted also as solutions of (H|) on a circle of length L = m5, for 
any m £ N. The integer m parameterizes different topological sectors: in the m-th 
the pendula chain twists around the circle m times. 

2 Explicit travelling-wave solutions with unit 
velocity 

The purpose of this work is to determine in closed form the travelling- wave solutions 
flU) just described. We first transform eq. (JSJ), with the help of the identities (fTBj) . 
into 

, F' A F(l-F 2 ) 

An = t — 4 — 

1+F 2 7 (1 + F 2 ) 2 

by looking for g in the form g = 4 tan -1 F and then into 

2ay' = 2y + 7 (l+y 2 ) (8) 

by looking for F in the form F = y+^Jl+y 2 . Note that diverging of \y\ at some point 
£o does not affect the continuity (and smoothness) of g at £o, even if the right limit 
is oo and the left one is — oo, or viceversa: y — * ±oo respectively implies F — > oo,0 
whence g mod 2tt in either case, which is compatible with a continuous g. 

Below we solve for ?/(£) explicitly. Putting all redefinitions together, we shall 
find solutions (p through the formula 



V9 ± (x,t)=4tan" 1 y(±x-t) + y/l+y 2 {±x-t) -vr. (9) 

Only if 7 < 1 the solutions y± = -^f~ x ± \J 7 -2 — 1 of the second degree equation 
y 2 +y2/7 + l = are real and therefore give (real) constant solutions = y± of 
([8]), whence the already mentioned constant solutions tp s ,(p u of ([1]). For nonconstant 
solutions (jHJ) is equivalent to 

= — — ri ( 10 ) 

7 y 2 + p + l 

separately in each interval where y' keeps its sign. The discussion of (fl~Qj) depends 
now on the value of the discriminant A = 4/7 2 — 4 of the equation y 2 +y2/7+l = 0. 
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If 7<1, then A > 0, y± are real and different and (|10p can be written as 



2a 



dy 



a 



7 {y-y+)(y-y-) y/i^f 



dy 



dy 



y-y+ y-y- 



which is integrated to give the two families of solutions 



2/1 (0 



y+ + y-e 



1 + e A($-to) ' 



2/2(0 



1 _ e A(« ) 



(11) 



where A := ( \A — 7 2 )o _1 and £o is an integration constant. One easily checks that 
y[,y 2 (and therefore also g[,g' 2 ) are respectively negative-, positive-definite; and 
that limg^-i-oo yi(£) = y=p for both i = 1,2. Using formulae (|23M24p shown in the 
Appendix we thus find 



lim F[ yi (£)] =F{yJ) = tan0, 

4— >oo 



.lim F[ yi (0}=F(y^ 



,'7T 

tan ( — - 
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for both i = 1,2, and mod 27r on one side a stricly decreasing g\ (£) with 



lim gi = 7r — sin 1 7, 

£^—00 



lim 51 = sin 1 7, 

£^00 



and on the other a stricly increasing 52(C) with 



lim Q2 = vr — sin 1 7, 



lim 52 = 27T + sin 1 7. 



As already noted, the singularity of 2/2 at £ = £ does not affect the continuity (and 
smoothness) of g 2 - Correspondingly, mod 2ir 



lim </3 1 
lim Lpt 



sin 1 7 = 93 s , 



■ — 1 s 
8111 7 = if 



lim <pf = -7r + sin _i 7 = 99", (12) 
x— »±oo 

lim (ft- = tt + sin -1 7 = </5 n , (13) 



therefore 99^, ^ are unstable solutions, as noted in [5]. 

If 7 = 1, then A = 0, y± = — 1 and (llOj) can be written as 



d£ = 2a 



(y + i) 



-2a d 



1 



y + 1 



which is integrated to give 



y(0 



This implies, with the help of 



lim = -1, 

£— >±oo 



1 + 



2a 



lim F[y(e)] = v / 2-l = tanJ, 



(14) 
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whereas again the singularity of y at £ = £o does not affect the continuity of g. As 
y', and therefore also F',g', are positive-definite, one finds mod 2ir 

lim g = — , lim g = — 

£— >— oo 2 £^oo 2 

and, correspondingly, 

_i_ 7T i 37T 

hm (p = --, lim = — ; (15) 

x— >:foo z x— >±oo z 

also these (/J 1 * 1 are unstable, as noted in [5]. 

Finally, if 7> 1, then A<0, y± are complex conjugate and the denominator of 
(fTUl) does not vanish for any value of y. Setting w := (y7+l)/\/7 2 — 1 ([10]) can be 
written as 

2a 



which is integrated to give £— £o = 2a tan 1 w/^^f 2 — 1, whence 

(16) 



y(0 = -7 1 + 1 = tan 



vT- 7 - 



2a 



where £o is an integration constant. This is a periodic function with period 

„ 27ra 



(17) 



and the latter is also the period occurring in ([6]). In fact, choosing £o = f° r 
simplicity, we see that as £ varies from —3/2 to E/2 varies from — oo to oo, 
F(£) varies from to oo, <?(£) varies from to 2tt. By continuity of g (which again 
is not affected by the singularity of y at £ = H(fc + l/2) (A; 6 Z), we thus find 
the behaviour ©. The corresponding solutions Cp^ fulfill ([7]), describe arrays of 
evenly-spaced kinks (see fig. [2]) moving with velocity ±1 and are stable [5] (see also 

13 El)- 

Appendix 

We first recall the trigonometric identities 

2tana 1— tan 2 a 4tana(l — tan 2 a) 

sin2a= 5—, cos2a= 5— sm4a = — ^—75 — . 

l+tan z a l+tan z a (l+tama) z 

(18) 

Given 7 e[0, 1], let 6 := \ sin _1 7 € [0, f ]. Then 7 = sin 40, \A-7 2 = cos 46* and, 
using the bisection formulae, 



l + \A-7 2 = V1+COS46* = cos 20, (19) 
1- \A-7 2 = V1-COS40 = \/2sin20, (20) 
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whence in turn 



V^-VI + aA-T 2 = V2(l-cos20) = 2\/2sin 2 6>, 



TT 



l-cos[ --20 



V2-\/l- \/l-7 2 = \/2(l-sin20) = \/2 
Hence, using also the sinus duplication formula, we end up with 



TT 

M,1_ ' 4" 



(21) 
(22) 





-7 2 


7 




2sin 2 (f- 




cos 26* 
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y/2-Jl-y/T- 



4 sin 20 sin 2 (f 



-0) 



sin 



sin 40 



tan 



cos 



TT 



72-vi+Vi-f 



4cos20sin 2 2sin 2 £ 



sin40 



sin20 



(23) 
tan0. (24) 



If we choose 7 = 1 in (|2ip and use the sinus duplication formula we find as 
another consequence 



V2-1 = 2\/2sin 2 
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